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We present the high time resolution observations of one 
of the Langmuir wave packets obtained in the source region 
of a solar type III radio burst. This wave packet satisfies 
the threshold condition of the supersonic modulational in- 
stability, as well as the criterion of a collapsing Langmuir 
soliton, i.e., the spatial scale derived from its peak inten- 
sity is less than that derived from its short time scale. The 
spectrum of this wave packet contains an intense spectral 
peak at local electron plasma frequency, / pe , and relatively 
weaker peaks at 2 f pe and 3 f pe . We apply the wavelet based 
bispectral analysis technique on this wave packet and com- 
pute the bicoherence between its spectral components. It 
is found that the bicoherence exhibits two peaks at (~ / pe , 
~ fpe ) and (~ f pe , ~ 2/pc), which strongly suggest that the 
spectral peak at 2/ pc probably corresponds to the second 
harmonic radio emission, generated as a result of the merg- 
ing of antiparallel propagating Langmuir waves trapped in 
the collapsing Langmuir soliton, and, the spectral peak at 
3 fpe probably corresponds to the third harmonic radio emis- 
sion, generated as a result of merging of a trapped Langmuir 
wave and a second harmonic electromagnetic wave. 


1. Introduction 

Solar flares, which represent the most dramatic energy 
releases from the Sun eject electrons into the corona and 
interplanetary medium. These electrons form bump-on-tail 
distributions and drive Langmuir waves unstable. These 
Langmuir waves subsequently couple to type III radio emis- 
sions at the fundamental and second harmonic of the elec- 
tron plasma frequency, /pe- This scenario, generally known 
as the plasma hypothesis [Ginzburg and Zheleznyakov , 1958] 
is supported by the observed fast negative frequency drift of 
the type III radio burst, which is interpreted as due to the 
electron beam propagating radially outward into lower and 
lower densities of the corona and interplanetary medium. 
Furthermore, the in situ observations of energetic electrons 
[Lin et al ., 1973, 1981; Ergun et al, 1998] and Langmuir 
waves [ Gurnett and Anderson , 1976, 1977; Lin et al ., 1986; 
Kellogg et al ., 1992; Gurnett et al , 1993; Thejappa et al. , 
1993; Hospodarsky and Gurnett , 1995; Thejappa and Mac- 
Dowall , 1998; Thejappa et al , 1999; Henri et al , 2009] in 
the source regions of type III radio bursts have provided 
additional support for the plasma hypothesis. Some times, 
faint third harmonic emissions are also observed in type III 
bursts [Benz, 1973; Takakura and Yousef, 1974]. One of 
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the outstanding questions is the identification of the mech- 
anism/mechanisms responsible for conversion of Langmuir 
waves into electromagnetic radiation. 

Some authors [Papadopoulos et al, 1974; Smith et al, 
1979; Goldstein et al, 1979; Goldman et al, 1980; Nicholson 
et al, 1978; Papadopoulos and Freund, 1978] predict that 
in type III radio bursts, the Langmuir waves are strongly 
turbulent, and therefore, the related oscillating two stream 
instability (OTSI) and Langmuir collapse [ Zakharov , 1972] 
play significant roles in the stabilization of electron beams as 
well as in conversion of Langmuir waves into electromagnetic 
waves. The Fast Envelope Sampler of Ulysses URAP exper- 
iment [Stone et al, 1992] have provided some evidence for 
strong turbulence processes in type III burst sources [Kel- 
logg et al, 1992; Thejappa et al, 1993, 1996; Thejappa and 
MacDowall, 1998; Thejappa et al, 1999; Thejappa and Mac- 
Dowall, 2004]. The in situ observations from the improved 
Time Domain Sampler (TDS) of the STEREO/WAVES ex- 
periment [Bougeret et al, 2008] (improved over that of all 
similar high time resolution receivers flown in earlier space- 
craft in precision, linearity, sample length and rate [Kel- 
logg et al, 2009]) and newly available higher order spectral 
analysis techniques have conclusively shown that the OTSI 
and related strong turbulence processes are commonly oc- 
curring phenomena in type Til burst sources [Thejappa et 
al, 2012a, b, c]. 

In this paper, we present new observations of a Langmuir 
wave packet captured by the STEREO TDS in the source 
region of a local type III burst. This wave packet satisfies 
the threshold condition of the supersonic modulational in- 
stability, as well as the criterion of a collapsing soliton; the 
spatial scale derived from its peak intensity is less than the 
spatial scale derived from the observed time scale. The spec- 
trum of this event, which is the focus of this study contains 
peaks at 2/ pe and 3/ pc in addition to an intense peak at / pe . 
To investigate whether any phase coherence exists between 
these spectral components, we compute the bicoherence us- 
ing the wavelet based bispectral analysis technique. This 
bicoherence exhibits two peaks, which strongly suggest that 
(1) the 2/ pe spectral peak probably corresponds to the elec- 
tromagnetic waves generated as a result of the coalescence 
of anti-parallel propagating Langmuir waves trapped in the 
collapsing Langmuir soliton, and (2) the 3 f pe peak corre- 
sponds to third harmonic radio emission generated as a re- 
sult of coalescence of trapped Langmuir wave with second 
harmonic electromagnetic wave. In section 2, we present the 
observations, in section 3, we present the bispectral analysis, 
and in section 4, we present the conclusions. 

2. Observations 

In Fig. 1, we present the dynamic spectrum of one of 
the local type III solar radio bursts (fast drifting emissions) 
and associated Langmuir waves (non-drifting emissions), ob- 
served on February 11, 2012 by the STEREO B spacecraft. 
In Fig. 2. we present the frequency- time spectrogram in a 
narrow frequency range, which shows the bursty structure 
of the Langmuir waves. The substantial frequency spread- 
ing of Langmuir waves seen here is probably due to non- 
linear frequency broadening. In this figure, the Langmuir 
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wave burst corresponding to the TDS event of current in- 
terest is shown by an arrow. In Fig. 3a, we present one 
of the most intense wave packets detected during this type 
III event. This wave packet is characterized by the peak 
electric field strength El of 64.6 mVm' 1 and ^-power du- 
ration D of ~ 11 ms. Since the E x and E z signals are 
weaker and show the same general features as the E y signal, 
we analyze only the E y signal. We justify one-dimensional 
treatment by assuming that these Langmuir wave fields are 
probably aligned along the ambient magnetic field. In Fig- 
ure 3b, we show the spectrum of this waveform from 12 to 
45 kHz. An intense peak at f pe ~ 14 kHz, and two sec- 
ondary peaks, one at ~ 28 kHz (second harmonic 2 fpe) 
and an another one at ~ 42 kHz (third harmonic 3/pe) > 
are clearly seen in this spectrum. The observed harmonics 
probably are the natural signals since: (1) the simulations 
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Figure 1 . Dynamic spectrum of the local type III ra- 
dio burst (fast drifting emission from ~ 3 MHz down to 
~ 20 kHz) and associated Langmuir waves (non-drifting 
emissions in the frequency interval 10-14 kHz.) 
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Figure 2. The frequency- time spectrogram of the Lang- 
muir waves observed during the type III burst of Febru- 
ary 11, 2012. The Langmuir wave emissions are very 
bursty and show substantial frequency spreading. The 
arrow shows the Langmuir wave burst corresponding to 
the current TDS event. 


of instrumental nonlinearities [Walker et a/., 2002] as well as 
the ground testing have shown that 3 f pe signal is stronger 
than 2 f pe signal [Malaspina et al, 2010], whereas in this 
case the higher the harmonic, weaker the signal, and (2) 
according to Kellogg et al. [2010], the STEREO/ WAVES 
experiment is designed in such a way that the responses in 
the Langmuir wave frequency range are as accurately lin- 
ear as possible, and no definite relation exists between the 
harmonic and fundamental amplitudes. 

Assuming that the frequency of main spectral peak in Fig. 
3b corresponds to the / pe , we estimate the electron density 
n c during this event as ~ 2.4 x 10 6 m -3 using the relation 
f P e = 9 n\‘ 2 . The solar wind speed v sw is measured as ~ 500 
kms -1 by the STEREO PLASTIC experiment [Galvin et al . , 
2008] , and we assume that the electron temperature T e dur- 
ing this event is ~ 10 5 K. We estimate the speed of the 
electron beam Vb as ~ 0.28c by fitting the frequency drift 
curve to the dynamic spectrum of the type III event (c is 
the speed of light). Using these values, we estimate: (1) De- 
bye length \De — 69T e 1/2 n7 1/2 ~ 14 m for T e ~ 10 5 K and 
n e ~ 2.4 x 10 6 m“ 3 , (2) wave number of Langmuir waves, 
~ 10 -3 m -1 for f pe ~ 14 kHz and Vb ~ 0.28c, 
(3) kL^De ~ 1.4 x 10 -2 , (4) normalized peak energy den- 
sity, — 5.6 x 10 -3 for El = 64.6 mVm -1 , 

rie = 2.4 x 10 6 m -3 and T e = 10 5 K, and (5) the j level 
spatial scale S ~ Dv $ w ~ 393A De for the measured values 
of D = 11 ms and v S w = 500 kms -1 . 
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Figure 3. (a): One of the most intense Langmuir wave 
packets captured by the Time Domain Sampler (TDS) 
during the type III event of Fig. 1. The ^ power du- 
ration of 11 ms which is equhalent to the spatial scale 
of 393 A^e is also shown, (b): The FFT spectrum of 
the TDS event, which shows the peaks at local electron 
plasma frequency, f pe as well as at 2 f pe and 3/ pe . 
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The threshold for the supersonic modulational instability 
[Zakharov, 1972] 
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is easily satisfied in this case, since ~ 5.6 x 10 3 , 

^ ~ 5.5 x 10 -4 and (k^De ) 2 ~ 2 x 10 -4 . The wave 
packet also satisfies the criterion of the collapsing envelope 
soliton [Thornhill and ter Haar , 1978; Gumett et al. , 1981] 
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where E denotes the expected value, X(f) is the Fourier 
transform of x(t) and * denotes complex conjugation. This 
equation clearly indicates that the bispectrum is zero unless 
three waves at frequencies fi , and fi -I- /2 are present in 
the time series with phase coherency amongst the waves. For 
spontaneously generated non-interacting modes, the phases 
would be random and therefore the statistical averaging will 
be zero. The bicoherence, which is the normalized bispec- 
trum, provides the quantitative measure of the phase coher- 
ence. It is defined as: 

}?(f f\ \BUhMt rn 

Ul,/J [E[\X{h)X{h)X-{h+hW W 


where Ak = ^ is the wavenumber characteristic of the en- 
velope, since the observed ~ 5.6 x 10 -3 is greater than 
(A/cAdc ) 2 ~ 2.6 x 10"" 4 estimated for the spatial scale of 
S ~ 393A De- 

3. Bispectral Analysis 

For nonlinear wave- wave interactions, the frequencies (/ 1 , 
/ 2 , fa) and wave vectors (fa, fe, fajjnust satisfy the reso- 
nance conditions: /i + / 2 = fa and /ci + fc 2 = • Since the 

wave number ki is related to the phase <j>i of the wave, the 
phases of three waves are also related as (pi 4- </> 2 = 03 • Bis- 
pectrum, a higher order spectrum, permits the distinction of 
the spontaneously excited normal modes from the coupled 
modes by measuring the degree of phase coherence between 
the modes. The FFT based bi-spectrum is defined as 

Btfuft) = E[X(fi)X{f 2 )X(fi+f2) m ]> (3) 
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Figure 4. (a): The wavelet based spectrogram of the 
TDS event of Fig. 3a, which shows the spectral features 
at f pe , 2/pe as well as at 3 f pe . (b): The wavelet based bi- 
coherence spectrum, which shows peaks at (~ 14, ~ 14) 
kHz and (~ 14, ~ 28) kHz. These bicoherence peaks sug- 
gest that their spectral components interact to produce 
the emissions at 28 kHz and 42 kHz, respectively. 


A unit value for the bicoherence indicates perfect coupling, 
a zero value indicates no coupling, and any value between 
zero and one indicates partial coupling. Periodogram esti- 
mation (interval segmentation) is usually used to estimate 
the bispectrum. The approximate variance of the bicoher- 
ence estimator is defined as [Kim and Powers , 1979]: 

var (b 2 )~^[l-b 2 Ui,f2)], (5) 

where N is the number of total data points and M is the num- 
ber of data records into which the time series is divided. To 
assess the low bicoherence levels, one needs long and station- 
ary time series with N >> 1. The bispectral analysis has 
been applied on the data from laboratory [Kim and Powers , 
1979] as well as space plasma experiments [Lagoutte et al. , 
1989; Dudok de Wit and KrasnoseVskikh , 1995; Bale et al. y 
1996; Balikhin et al. , 2001; Walker et al. , 2003; Henri et 
al , 2009]. 

The wavelet based bispectral analysis technique is su- 
perior to the FFT based technique, since it can handle 
both short and non-stationary time series. It has been 
developed to study the strong turbulence processes at the 
Earth’s quasi-parallel bow shock [Dudok de Wit and Kras- 
noseVskikh , 1995], and it has been successfully used to study 
the electrostatic decay in type III sources [Henri et al., 2009]. 
In this study, we apply this wavelet based bispectral analy- 
sis technique on the wave packet of Fig. 3a. The continuous 
wavelet transform (CWT), which is used in this technique 
is defined as 

W{a,r) = -J= 

vM 

where the scale a is related the frequency / as a = j and 
r is the time. The Morlet wavelet xp(t) = e l2,r V 2 is gener- 
ally used to study nonlinear wave phenomena. The wavelet 
based bispectrum in the frequency space can be defined as 

B(/i,/ 2 ) = E[W(f l )W(h)W*(f l +/ 2 )], (7) 

where W(n8t,f = 1/a) is a function of time t as well as 
frequency f. The wavelet bicoherence can be defined as 

l2/ r r\_ \B(f\,f2)\ 2 /«x 

Ul,/2; (E[\W(fi)W(f 2 )W*(f 1 +f 2 )\]) 2 ' W 

The discrete wavelet coefficients W x of the signal are com- 
puted using the Morlet wavelet. The discretization frame 
of the wavelet coefficients is taken regular with a constant 
frequency step to ensure the frequency resonance conditions. 

In Figure 4a, we present the modulus of the wavelet trans- 
form applied to the waveform data presented in Fig. 3a. The 
spectral features corresponding to the fundamental, second 
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and third harmonic components are clearly seen in this spec- 
trogram, which show that the higher the harmonic, weaker 
the emission. It also shows the temporal coincidence of the 
second and third harmonic emissions with collapsing Lang- 
muir soliton. In Figure 4b, we present the results of wavelet 
based bicoherence analysis. This bicoherence spectrogram 
clearly shows evidence for three wave- wave interactions. The 
bicoherence peaks at (~ 14 kHz, ~ 14 kHz), and (~ 14 kHz, 
~ 28 kHz) suggest that these spectral components inter- 
act to produce the emissions at 28 kHz and 42 kHz, re- 
spectively. The most probably three wave interactions are 
L + L -► T 2fpe , and L + T 2fpe -> T 3/pe , respectively (L and 

L are the opposite propagating Langmuir waves, and T 2 / pc 
and T 3 / pe are the second and third harmonic electromag- 
netic waves, respectively). 

The Langmuir waves involved in the interaction L + 
L — * T 2/pe are probably the anti-parallel propagating waves 
trapped in the collapsing soliton. This is consistent with the 
wavelet spectrogram (Fig. 4a). which shows the temporal 
coincidence of the second and third harmonic emissions with 
the peak of the Langmuir wave packet. The wave- wave in- 
teraction L+L — ♦ T 2 / pc agrees with that of Papadopoulos et 
al. [1974], who suggested that the anti-parallel propagating 
modes generated during oscillating two stream instability 
(OTSI), get trapped in the collapsing soliton and couple to 
each other generating electromagnetic waves at 2/ pe . The 
interpretation of the spectral peak at 28 kHz in terms of 
T 2 / pe is consistent with that of Bale et al. [1996], who in- 
terpreted the harmonic spectral peak of the wave packet 
observed by the WIND TDS in the upstream solar wind in 
terms of T 2 / pc . One should note that the observed phase 
coherency between the fundamental and second harmonic 
spectral components does not support the mechanism in- 
voked by Papadopoulos and Freund [1978] as well as Gold- 
man et al. [1980], which invoke the spontaneous emission of 
electromagnetic waves at 2 f pe by the stable and collapsing 
solitons, respectively; the primary waves at f pe and harmon- 
ics at 2 fpe in such a case are not phase correlated. 

Malaspina et al. [2010] have interpreted the harmonics 
of the wave packets captured by the STEREO TDS in 
the solar wind in terms of localized nonlinear electrostatic 
modes driven by eigenmode-localized Langmuir waves. Fur- 
thermore, these authors have suggested that these nonlin- 
ear electrostatic modes probably emit the coherent electro- 
magnetic radiation spontaneously at respective frequencies. 
However, the observed phase coherency between the funda- 
mental and harmonic components rules out (1) the possi- 
bility that the 2 f pe spectral peak corresponds to harmonic 
Langmuir waves, since the non-linear wave- wave interaction 
between the fundamental and harmonic Langmuir waves 
is forbidden (harmonic Langmuir waves are not the eigen 
modes), and (2) the second harmonic radio emission is not 
due to any antenna emission but is due to wave-wave inter- 
actions between the opposite propagating Langmuir waves. 

Kellogg et al. [2010] argued that the harmonics of the 
wave packets captured by the STEREO TDS in the Earth’s 
fore shock are probably due to the electrons trapped in the 
Langmuir wave potential maxima. In this scenario also, the 
bicoherence which measures the phase coherency is expected 
to be zero contrary to the observed values, which suggests 
that the oberved harmonic modes in the present case are not 
due to trapping of electrons in the Langmuir wave maxima. 

As seen from Figs. 3b and 4a, the third harmonic emis- 
sion is usually very weak, and, therefore, there are only a 
handful of observations of third harmonics in solar radio 
bursts. Zheleznyakov and Zlotnik [1974] have suggested two 
kinds of wave- wave interactions, namely, L + L + L —> Tzf pe 
or L + T 2 f pe —* T 3 f pe as possible mechanisms for excitation 
of the third harmonic. However, the observed bicoherence 
peak at (/ pe , 2 f pt ) strongly suggests that the third harmonic 
is probably generated by the later mechanism. 


4. Conclusions 

We have presented the high time resolution observations 
of a Langmuir wave packet associated with a local solar type 
III radio burst. This wave packet is found to satisfy not only 
the threshold condition of the supersonic modulational in- 
stability, but also the criterion for a collapsing Langmuir 
soliton. Furthermore, the spectrum of this wave packet is 
found to contain an intense peak at f pe and weaker peaks at 
2 f pe and 3/ pe . Since the fundamental is more intense than 
the second harmonic which is stronger than the third har- 
monic, contrary to the results obtained by the ground testing 
and simulations of instrumental effects which have predicted 
a stronger third harmonic in comparison to the second, it is 
argued that the observed harmonic peaks probably are not 
due to instrumental effects. In order to test whether the 
fundamental and harmonic spectral components are phase 
coupled with each other, the wavelet based bicoherence anal- 
ysis technique is applied on this wave packet. It is shown 
that the bicoherence exhibits two peaks at (~ f pe = 14 kHz, 
~ fpe = 14 kHz) and (~ f pe = 14 kHz, - 2 f pc = 28 kHz). 
The bicoherence peak at (~ 14 kHz, ~ 14 kHz) strongly 
suggests that (1) the three wave- wave interaction ~ 14 kHz 
+ ~ 14 kHz — 28 kHz, i.e., L + L — > T 2 / pe occurs in the 
regime of strong Langmuir turbulence, (2) the anti-parallel 
propagating Langmuir waves involved in these wave-wave 
interactions are probably generated by the oscillating two 
stream instability (OTSI), which is responsible for the for- 
mation of collapsing Langmuir solitons as suggested by Pa- 
padopoulos et al. [1974], (3) the observed second harmonic 
spectral component is not due to the spontaneously emitted 
coherent emission by the stable and collapsing solitons as 
suggested by Papadopoulos and Freund [1978] and Goldman 
et al. [1980], respectively, because in the case of coherent 
emission, the bicoherence is expected to be very close to 
zero, and (4) the harmonic modes are not the electrostatic 
modes spontaneously generated by either the eigenmode- 
localized Langmuir waves as suggested by Malaspina et al 
[2010], or by the trapped electrons in the Langmuir wave 
potential maxima as suggested by Kellogg et al. [2010]. As 
far as the bicoherence peak at (~ f pe = 14 kHz, ~ 2/ pc = 28 
kHz) is concerned, it provides evidence for the merging of 
the Langmuir waves with the second harmonic electromag- 
netic wave yielding a third harmonic electromagnetic wave, 
i.e., L + T 2 f pc — ► T 3 f pe as suggested by Zheleznyakov and 
Zlotnik [1974]. 
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